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1.  Introduction.  There  are  at  least  five  different  formulations  of 
duality  --  the  original  Fenchel  formulation  [3,8],  the  (generalized) 
geometric  programming  formulation  [1,U,5],  the  Fenchel-Rockaf ellar 
formulation  [6,8],  the  ordinary  Lagrangian  formulation  [10,2,8],  and 
the  Rockafellar  formulation  [7,8,9].  Although  each  formulation  has 
its  own  advantages  and  disadvantages,  each  can  also  be  viewed  as  a 
special  case  of  each  of  the  other  four. 

The  appropriate  specializations  have  already  been  carried  out  in 
[ U ,8 , 9] , but  only  for  a very  limited  geometric  programming  formulation. 
Thus,  this  paper  complements  those  references  by  specializing  the  Rocka- 
fellar formulation  to  the  most  general  geometric  programming  formula- 
tion. , 

Section  2 presents  a version  of  the  Rockafellar  formulation  that 
facilitates  the  specialization  given  in  Section  3.  This  specializa- 
tion does  not  require  convexity  assumptions  and  uses  only  elementary 
real  analysis. 

2.  Rockafellar  duality.  Suppose  that  g:C  is  a (proper)  function  g 
with  a nonempty  (effective)  domain  C c , and  assume  that  the  inde- 
pendent variable  (d,p)  in  C is  the  Cartesian  product  of  a "decision" 
(vector)  variable  d and  a "perturbation"  (vector)  parameter  p . 

Consider  the  parameterized  family  G that  consists  of  the  follow- 
ing optimization  problems  A(p)  . 


I 


i 


1 


2 


. 


* 


1 


PROBLEM  A (p ) . Using  the  "feasible  solution"  set 

s(p)  = {d|(d,p)  <=  C]  , 

calculate  both  the  "problem  Infimum" 

<p(p)  = inf  g(d,p) 

d es(p) 

and  the  "optimal  solution"  set 

S*(p)  = {d  € S(p)|g(d,p)  = cp(p )}  • 

For  a given  perturbation  p , problem  A(p)  is  either  "consistent" 
or  "inconsistent",  depending  on  whether  the  feasible  solution  set  S(p) 
is  nonempty  or  empty.  The  (effective)  domain  of  the  infimum  function  cp 
is  the  "feasible  perturbation"  set 

P = {p|problem  A(p)  is  consistent]  , 

which  is  obviously  identical  to  [p|(d,p)  € C for  at  least  one  d]  and 
hence  is  not  empty.  Unlike  the  function  g , the  function  cp  may 
assume  the  value  -®  . However,  for  our  purposes,  it  is  not  advan- 
tageous to  follow  Rockafellar ' s custom  of  artificially  defining  g and 
Cp  to  be  +oo  outside  their  respective  domains  C and  p . 
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Now,  suppose  that  g:C  has  a "conjugate  transform"  h:D  ; that  is, 
suppose  there  is  a function  h with  a nonempty  domain 

D = {(q,e)|  sup  [<q,d)  + <e,p)  - g(d,p)]  < +®} 

(d,p)€C 

and  function  values 

h(q,e)  = sup  [<q,d)  + <e,p)  - g(d,p)]  . 

(d,p)€C 

The  inner  product  (q,d)  associates  the  "dual  perturbation"  parameter  q 
with  the  "primal  decision"  variable  d , and  the  inner  product  (e,p) 
associates  the  "dual  decision"  variable  e with  the  "primal  perturbation" 
parameter  p . 

Consider  the  parameterized  family  KJ  that  consists  of  the  following 
optimization  problems  B(q)  . 

I 

PROBLEM  B ( q ) . Using  the  feasible  solution  set 

T(q)  = {e|(q,e)  £ D)  , 
calculate  both  the  problem  infimum 


\(r(q)  = inf  h(q,e) 
« €T(q) 


- I*  - 


and  the  optimal  solution  set 


T*(q)  = {e  € T(q)|h(q,e)  = \|r(q)}  . 

Needless  to  say,  the  domain  of  the  infimum  function  is  the 

feasible  perturbation  set 

Q = {q|problem  B(q)  is  consistent]  , 
which  is  clearly  not  empty. 

Due  to  the  known  symmetry  [3,8]  of  the  conjugate  transformation 
on  the  class  of  all  closed  convex  functions  g:C  (as  well  as  the 
obvious  symmetry  of  the  preceding  association  of  perturbation  parame- 
ters and  decision  variables),  families  G and  B are  termed  Rocka- 
fellar  dual  families,  and  problems  A (0)  and  B(0)  are  termed  Rocka- 
fellar  dual  problems.  Actually,  Rockafellar  [7,8,9]  formulates  B 
as  a family  of  maximization  problems  by  placing  minus  signs  in  front 
of  the  sup  and  e in  the  definition  of  h:D  . Although  that  formula- 
tion facilitates  specializations  [7,8]  to  (the  standard  formulations  of) 
linear  programming  duality  and  ordinary  programming  duality,  the  pre- 
ceding formulation  will  facilitate  our  specialization  to  geometric  pro- 
gramming duality. 

To  (re)orient  the  reader  toward  the  preceding  formulation,  we  now 
summarize  Rockafellar ' s main  results  in  terms  of  that  formulation.  In 
particular,  the  primal  Infimum  function  q>  is  finite  everywhere  on  its 


domain  P and  has  a con  jugate  transform  (with  a nonempty  domain  and 


and  suppose  that 


C = { (x,x,u,n)|xk+uk  € c^,  k€{o)UI ; (x^+u^ ,k ^ )£ c*,  J € J ; 


x6X;  and  g^(x  +u  ) + ^*0,  i 61} 


/ > A _ / » A / 0 0 \ + / j J \ 

g(x,K,u,n)  = G(x+u,k)  = gQ(x  +u  ) + 2*  8j(x  +uJ,Xj) 


Then,  the  Rockafellar  family  G is  the  "geometric  programming 
family"  F (described  in  subsection  3.3.5  of  [5]);  and  the  crucial 
question  now  is  whether  the  Rockafellar  dual  family  fi  is  the  "geo- 
metric programming  dual  family"  G (described  in  subsections  3.3.4  and 
3.3.5  of  [5]).  To  obtain  the  answer,  we  need  to  compute  the  conjugate 
transform  h:D  of  g:C  in  terms  of  both  the  "dual"  Y of  the  given 
cone  X and  the  conjugate  transforms  h^rD^  °f  the  given  functions 

gk:Vk€f(y>UlUJ  ’ 

To  compute  h:D  , assume  that 


q = (v,v)  and  e = (z,\)  , 


where  v has  the  same  component  partitioning  as  x , and  where  z has 
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I 


the  same  component  partitioning  as  u . Then 


h(v,v,z,\)  = 8uP  { (v°,x°>  +V  (v1^1)  + 53  <v^x^>  + £ vk. 

(x,K,u,n)  T J J J J 

< 0 0.  / i i.  / j j . «— i / 0 0 > 

+ <z  ,u  > + £ <z  ,u  ) + E<z  >u  > + E Xj 4. "8n(x  +u  ) 

I J I 

- E 8j  (xJ+uj  ,k^  ) |x°+u°  € CqJx^u1  € C±,i  g I ; (xJ+uJ  ,Kj  )g  C* , j £ J 


x € X;  and  g^x^+u1)  + ^ £ 0,  if  I], 


which  is  clearly  finite  only  if  X ;>  0,  i £ I ; in  which  case 

h(v,v,z,X)  = sup  { [<v°,x°) + <z°,u°)  - g (x°+u0)] 

(x,/c,u)  0 

+ E C^vi<xi)  + - X g (x^u1)] 

I 1 1 

+ E Kv^,xl>  + (z^,u^>  +VjKj-gj(x^+u-',K:j)]  | 

x°+u°  ec0;xi+ui  €Ct,i  ei;(xJ+uj,^)  €C+,  j £ J;x  £X] 
or 


h(v,v,z,X)  = sup  { [<z0+v°,x0+u0>  - g (x°+u°)] 

(x,K,u)  0 

E,i  i i i i.  . * i i , . 

[(z  +v  ,x  +u  ) - X.g, (x  +u  )] 

I 1 1 

+ El^+v^x^+u-5)  +vj'<j  " gj(xJ+uj,Afj  )] 

•J 

- <z,x>  - <v,u)|x°+u0  e CQ ;x1+u1 6 C1,i g I; 


(xJ+uJ,<.)  ect,j  € J;x€X]  . 


The  following  lemma  provides  another  condition  that  is  necessary 
for  the  finiteness  of  the  preceding  expression. 


Lemma  A . The  preceding  expression  for  h(v,v,z,X)  is  finite  only  if 
z - v $ Y . 


Proof,  if  z-v^Y  , then  there  is  an  x fX  such  that  (z-v,x)  < 0 ; 
in  which  event  we  choose  c so  that:  c £C^;  c i £ I ; and 

(c^Aj)€Cj>  j€J,  f°r  some  fixed  x s 0 . Letting  x(s)  = sx  and 
u(s)  = c - sx  , we  observe  that  (x( s ) ,k , u( s ) ) satisfies  the  restric- 
tions on  (x,(c, u)  for  each  s 2 0 . Thus 


r r/  u u -O.  ,-0,.  V'  r/  1 i -iv  , 

{[<z  +v  ,c  > - gQ(c  )]  + 2^[<z  +v  ,c  ) - Xigi(c  )] 


+ £[<ZJ+VJ,CJ>  + VjKj  - gj(cj,Kj)]  - <Z,SX>  - <V,C-SxS} 


s h(v ,v,z,X) 


and  hence  h(v,v,z,X)  = +0°  because 


lim  {-(z,sx)  - (v,c-sx)}  = lim  f-(v,c)  + s(v-z,x)}  = +® 


by  virtue  of  the  property  (z-v,x^  < 0 


q . e .d  , 


Now,  if  z-v£Y  , then  z = y+v  for  some  y 6 Y ; in  which  case 
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h(v,v,z,X)  = 


rr/0  0 0 0«  , 0 0Nl 

sup  {[<z  +v  ,x  +u  ) - g (x  +u  )] 

(x,<,u)  0 

Lr/  i i 1 i,  , , i in 

[(z  +v  ,x  +u  > - X.g. (x  +u  )] 

I 1 


+ 5J[<zj+vJ,xj+uJ>  + VjKj  - gj(xJ+UJ,Xj)] 


• (y>x)  ' <v,x+u> |x°+u° £ CQ;  xi^ui(;C1,  ifl  ; 


(x^+u^  ) € Cj  , J€J  ; XU, 


Since  0 £ <y,x>  for  each  xfX  , it  is  easily  seen  that 


h(v,v,z,X)  = 


sup  {[<z°+v°,u0>  - g (u0)] 

(<,u)  0 

+ 13^zi+vi--ui^  * Xigi(u1)] 


+ £ [<ZJ+Vj,U^  + Vj(C  j - gj(uJ,Kj)]  - <v,u)  | 

u°<=C0;  ui£Ci,  i € I;  (uJ,Kj)€C+J  J € J} 
= { nSuP  [(z°,u°>  - g0(u°)]  + £ .sup  [(z^U1)  - x.g.  (u1)] 

..U  c n u t r n 


uufC 


I u £ C 


+ £ <sup  . KzJ,uJ>  + v.k.  - g^(uJ,«  ; 
J J J j J 


Consequently,  (v,v,z,X)gD  if  and  only  if:  X^O,  i € I S z = y+v  for 

some  y $ Y ; and  each  term  on  the  right-hand  side  of  the  preceding  equations 
is  finite.  Of  course,  the  first  term  is  finite  if  and  only  if  z°  £ DQ  , in 
which  case  the  first  term  is  equal  to  h^(z^)  . The  finiteness  of  the 
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remaining  terms  can  be  conveniently  characterized  with  two  lemmas. 

The  following  lemma  characterizes  the  finitenes6  of  the  terms 
involving  the  index  set  I. 


Lemma  B 


. Given  that  \±  2 0 , the  |up  [(z^uS  - X^^u1)]  is  finite  if 


C 


^ i 

and  only  if  (z  jX^D*  , in  which  case 


sup 
u1  gc 


[(z1^1)  - Xtgi(u1)]  = h^(zt,Xi)  • 


Proof.  Simply  observe  that 


sup  (z^u1)  if  \t  = 0 


u1  6C, 


sup  [(zV1)  - X^Cu1)]  = | X^^zVx^  if  Xt  > 0 and  z1  6 X^ 


ul  gc 


+ °° 


if  Xt  > 0 and  z*  £ XjD^, 


+ .r,+ 


and  then  use  the  defining  formula  for 


q .e .d . 


The  next  lemma  characterizes  the  finiteness  of  the  terms  involving 
the  index  set  J . 


Lemma  C.  The  sup  [(z^,u^)  + “ g*(u^>K,)]  is  finite  if  and 

-JgD.  and  hj(zJ)  + v.  * 0 , in  which  cas_e 


only  if  both  zJ£Dj 


sup  [<ZJ,UJ)  + VjKj  - gj(u^,<j)l  - 0 . 
(uJ,icJ)  € 
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Proof.  First,  observe  that 


Sup  [(Z^U-*)  + Vj/Cj  ~ gj(^,Kj)] 

(u^x^eCj 

= sup  [sup{(z^,u^)  + Vj/Cj  - gj(u^ ) I (u^ >‘<j  ) € Cp  ] 
Kji 0 u^ 


= sup  [VjKj  + sup{<zJ,J>  - gj(uJ,Kj)|  (uJ,Kj)  gcp] 

■1^ 

/sup{<z^ ,u^)  - sup  (u^,d^>|  8 up  (u^,d^>  < +®} 

I -j  dJ0D, 


K^O 


= sup 
ICjSsO 


if  Kj  = 0 


viKi  + 


dJ^D 


\sup((zJ,uJ>  - Xjgj (uJ/k:j ) | uJ/#Cj  €Cj] 


.J  \ I .J 


If  1C  > 0 


= 8 up 

KyiO 


VJKJ  + 


/ 0 

+ 03 
+ «r 


if  Kj  = 0 and  z^  f D j ^ 

if  = 0 and  z^  £ 

if  Kj  > 0 and  z^  ft  Dj 


yKjhj(z-*)  if  Kj  > 0 and  z^D ‘ 


where  the  final  step  makes  use  of  the  fact  that  the  zero  function  with 
domain  D ^ (the  topological  closure  of  D ^ ) is  the  conjugate  trans- 
form of  the  conjugate  transform  of  the  zero  function  with  domain  . 

Now,  note  that  the  last  expression  is  finite  only  if  z^  € , in 

which  case  the  last  expression  clearly 


Sup  [VjKj  + K-jhj(zJ)]  . 


Kj20 
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But  this  expression  is  obviously  finite  if  and  only  if 
in  which  case  this  expression  is  clearly  zero. 


h («J)  + Vj  S 0, 

q.e.d. 


We  have  now  shown  that 

D = { | X±  i 0,  i €1 ; z=y+v  for  some  y€Y;  z°  €DQ; 

(zi,Xi)  €D^,i€I;  zj  €D.  and  hj  (zJ  Hv^O,  j 6J] 

and 

h(v,v^A)  = h0(z°)  + ]Ch^(zi,Xi)  • 

Consequently,  the  Rockafellar  dual  problem  B(0)  is  the  "geometric 
programming  dual  problem"  B (described  in  subsection  °f  [5])* 

Although  the  Rockafellar  dual  family  Hi  is  slightly  different  from  the 
geometric  programming  dual  family  G (alluded  to  in  subsection  3*3*5 
of  [5]),  the  difference  is  inconsequential  in  view  of  the  final  para- 
graph of  section  2 . 

However,  the  relation  z = y+v  shows  that  y can  be  used  instead 
of  z as  a dual  decision  variable  --  a change  of  variables  that  clearly 
induces  a one-to-one  mapping  from  B onto  G . In  particular,  this 
mapping  simply  translates  the  (domain  of  the)  dual  objective  function 
h(v,v;‘)  through  (*v,0)  --  a mapping  that  clearly  leaves  the  problem 
inf imum  \|r(v,v)  invariant  while  translating  the  optimal  solution  set 
T (v,v)  through  (-v,0)  . 
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